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We analyze theoretically and experimentally the far-held viscous decay of a two- 
dimensional inertial wave beam emitted by a harmonic line source in a rotating 
huid. By identifying the relevant conserved quantities along the wave beam, we 
show how the beam structure and decay exponent are governed by the multipole 
order of the source. Two wavemakers are considered experimentally, a pulsating 
and an oscillating cylinder, aiming to produce a monopole and a dipole source, 
respectively. The relevant conserved quantity which discriminates between these 
two sources is the instantaneous howrate along the wave beam, which is non-zero 
for the monopole and zero for the dipole. For each source the beam structure and 
decay exponent, measured using particle image velocimetry, are in good agreement 
with the predictions. 


I. INTRODUCTION 

In rotating fluids, the restoring action of the Coriolis force allows for the propagation of 
anisotropic, transverse, circularly polarized waves called inertial wavesi These waves are of 
fundamental interest for geo- and astrophysical flowss^^— they can for instance be excited in 
the fluid core of planets by tidal motions, precession or libration^^— Internal gravity waves 
in stratified fluids, relevant to the ocean and the atmosphere, share a number of properties 
with inertial wavesi^ Internal waves, or mixed internal-inertial waves when rotation and 
stratification effects are of comparable magnitude can also be excited in the ocean by 
the interaction of tides with topographyji^^— 

We focus in this paper on the scaling of the viscous decay of a wave beam emitted by a line 
source in a uniformly rotating fluid. Although this problem has received much attention, the 
influence of the multipole order of the source has not been addressed so far. Interestingly, 
whereas the growth of the wave beam thicknesses 5{x) ^ (with x the distance from the 
source), is independent of the multipole order of the source, the decay of the wave amplitude 
depends on which quantity is conserved along the wave beam (flowrate, momentum, or 
higher order moment), which is directly governed by the multipole order of the source. 
Such dependence was discussed in the case of internal waves produced by a point source 
in stratified fluids by VoisiniSS Here we address this problem for a two-dimensional inertial 
wave beam produced by a line source, extending the simple far-held quasi-parallel approach 
of Cortet et to a source of arbitrary multipole order. We show that the wave decay 
is steeper as the multipole order of the source increases, which implies that the far-held 
decay of a wave beam emitted by an arbitrary source is dominated by its lowest multipole 
component. 

Inertial wave beams emitted from localized sources are relevant to a broad range of labo¬ 
ratory and natural hows, including local forcing by a wavemaker immersed in the huid, but 
also global forcing such as precession, libration or tidal motion acting at the scale of the 
huid domain. This is because in all cases wave beams are emitted from critical lines, where 
the local slope of the solid boundaries equals the propagation angle of the wave. Along 
such critical lines the oscillating boundary layer erupts and forms oscillating beams in the 
bulk of the howj^ i^dd^^ .^ In conhned huid domains such beams rehect and, in the presence 
of sloping boundaries, may focus on wave attractorsThe eruption at critical lines 
produces two types of wave beams, associated with different magnitudes and scaling laws, 
propagating in planes tangent and non-tangent to the solid boundaryi^ Wave emission in 
the tangent plane, developing only along convex boundaries, is stronger: this is the case for 
conical inertial wave beams emitted from the inner core of a rotating spherical shellj^d^ or 
for internal wave beams excited by oceanic internal tides on the edge of continental shelves 
or ridges»i^“— Wave emission in the non-tangential plane, both at critical lines of convex 
or concave slope (for instance on the outer sphere of a rotating spherical shell^iid^ or at the 
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(a) pulsating cylinder 


(b) oscillating cylinder 




FIG. 1. (Color online) Sketch of the far-held wave beams resulting from a combination of parent 
beams emitted from the critical lines, for (a) a pulsating cylinder and (b) an oscillating cylinder. 
The size of the disturbance is assumed much larger than the thickness of the Ekman boundary layer. 
The dashed circles represent the separation between bimodal and unimodal regions. Eight wave 
beams are emitted from four critical lines (red points), which behave as local sources of non-zero 
instantaneous flowrate. Far from the cylinder, two parent beams propagating in the same direction 
combine and form a unique beam. The instantaneous howrate of the resulting merged beam is 
non-zero when the two parent beams are in phase (case a), and is zero when they are out of phase 
(case b). 


sloping bottom of a ridgei^), is of weaker amplitude. Emission of inertial wave beams is also 
observed from horizontal edges in containers such as a cylinder— >21 or a parallelepiped— 

A key property of an inertial wave beam spawned from an erupting boundary layer is its 
non-zero instantaneous flowrate: it can be modeled in the far field as originating from a 
monopole line source. Depending on the topology of the fluid domain, in particular on the 
distribution and relative phases between such elementary monopole sources, different wave 
beams propagating in the same direction may combine and form far-held beams of either 
non-zero or zero instantaneous flowrate, therefore corresponding to an effective monopole or 
a higher order multipole source. Note that such combination of beams requires propagation 
over a distance much larger than the separation between the sources, a requirement which 
is usually not satished in geo- and astrophysical situations. 

We restrict in the following to wave beams produced by effective line sources surrounded 
by the fluid. In practical situations, such line source corresponds to a two-dimensional convex 
disturbance, say a cylinder, dehning four critical lines (Fig. I). If the cylinder is pulsating 
(Fig. la), the periodic emission and suction of mass from these critical lines is in phase, so 
the far-held merged beams have non-zero howrate: this dehnes an effective monopole source. 
On the other hand, if the cylinder is oscillating (Fig. lb), the parent beams propagating 
in the same direction are out of phase, resulting in far-held merged beams of zero howrate: 
this dehnes an ehective dipole source. The pulsating and oscillating cylinders are therefore 
generic conhgurations to investigate the inhuence of the multipole order of a line source on 
the properties of the far-held wave beams. 

Most of the experiments investigating wave beams emitted by a local disturbance, both 
in stratihe d^^i^^~ — and rotatingiSiM huids, are based on oscillating wavemakers (with the 
exception of Makarov et al^ who report results from a pulsating cylinder in a stratihed 
huid). The resulting far-held wave beams are therefore distinct from those produced from 
an erupting boundary layer in a globally forced huid domain. The aim of the present paper 
is to compare the far-held properties of inertial wave beams emitted from a pulsating and an 
oscillating cylinder in a rotating tank, aiming to produce a monopole and a dipole source. 
Velocity measurements in the wave beam are achieved using two-camera multi-resolution 
particle image velocimetry, ensuring a good resolution both in the near and far helds. The 
two sources produce distinct decay exponents and wave beam prohles, in good agreement 
with the theoretical predictions. 
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II. THEORETICAL BACKGROUND 


A. Dispersion relation and viscous spreading 


The geometrical properties of inertial waves follow from their dispersion relationji 


cr = ± 


2rjk 

~W 


2ncos0, 


( 1 ) 


with cr > 0 the wave frequency and 6 the polar angle between the wave vector k and the 
vertical rotation vector Jl. Waves emitted from a localized harmonic disturbance propagate 
energy in directions making an angle ±0 to the horizontal, along two cones for a point source 
and along four plane beams for a line source normal to We restrict in the following to 
the line source configuration, for which the spatial decay of the wave is purely governed by 
viscosity. In each wave beam, fluid particles describe anticyclonic circular translations in 
the tilted plane normal to k. Since only the orientation of k is prescribed by the dispersion 
relation o but not its magnitude, the characteristic sizes of the wave (wave length, beam 
thickness) are governed by the boundary conditions and viscosity. 

The viscous spreading of the wave beam results from the combination of the energy 
propagation in the longitudinal direction x and its diffusion in the lateral direction z (see 
Fig. [2|). Its scaling can be obtained from a classical boundary-layer argument: During a 
time t, the wave energy spreads laterally over a distance (5 ~ with v is the kinematic 
viscosity, and propagates over a distance x = Cgt, where Cg is the group velocity. Evaluating 
Cg = {a/k)ta,n9 for the dominant wave number at a distance x from the source, k ^ 
simply yields 

6{x) - (2) 

where we introduce the viscous scale 

i = {vjatanOy^'^. (3) 


Although the scaling of the wave amplitude strongly depends on the multipole order of the 
source, the scaling of the beam thickness ([2]) is independent of the nature of the source, 
provided that 6{x) is much larger than the size of the source. 


B. Boundary layer equations 

We derive now the similarity solutions for a viscous 2D inertial wave beam emitted by a 
harmonic line source, focusing on the spatial decay of the wave amplitude and its dependence 
on the multipole order of the source. The derivation follows that of Thomas and Stevenson^ 
for internal waves in stratified fluids. We use the velocity-vorticity formulation of Cortet et 
al.^ generalized here to a source of arbitrary order. 

We consider a line source along the Y axis, of angular frequency cr, in a fluid rotating at 
rate D about the Z axis (Fig.[2|). The four wave beams emitted by the source being invariant 
along y, energy propagates in the (X, Z) plane. In the following, we consider only the wave 
beam propagating in one given quadrant. We start from the linearized vorticity equation in 
the rotating frame 


= (20-V)u-f (4) 

with u the velocity and a; = V x u the vorticity. We project (|H) on the local frame of 
the far-field wave beam, (e^,, ej^, e^), with aligned with the group velocity, making angle 
0 = cos“^(cr/2D) with the horizontal, and assume that the flow inside the wave beam is quasi¬ 
parallel (boundary layer approximation), i.e., such that |ua;|,|uy| |mz|; |a;a;|,|a;y| ^ 

and ~ &l. We introduce the complex velocity and vorticity fields in the tilted plane 
{x,y), U = Ux + iuy and W = uJx + itOy — idzU. Searching for harmonic solutions of the 
form U = Eq. (jH) writes 


dxUo -I- il^dlUo = 0. 


( 5 ) 
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FIG. 2. Two-dimensional wave beam emitted from a line source at {X = 0, Z = 0) in a fluid 
rotating about the Z axis. Only the beam propagating in the first quadrant {X > 0, Z > 0) is 
shown. The contour C shows the frontiers of the domain of integration D used in Sec. Ill Cl 


This equation admits similarity solutions as a function of the reduced transverse coordinate 
T] = of the form 


Uo{x,z) = Uo(^ fiv), (6) 

where Uq is a complex velocity scale and a > 0 the decay exponent to be determined. The 
solution considered in Cortet et a/.— was derived for the particular case a = 1/3.— Inserting 
Eq. dni) in Eq. ([5|) yields 


Sf" + ivf + 3iaf = 0. (7) 

Solutions of this equation, first given by Moore and Saffman^l for the problem of a vertical 
steady shear layer {9 = 7r/2), and later by Thomas and Stevensort^S in the context of internal 
waves, are 


pOC 

fmiv) = Cmiv) + (8) 

Jo 

with Cm and Sm real functions, even and odd respectively. The properties of these functions 
were considered in detail by Voisiniii Integrating Eq. m by parts gives 3fm+3 - ivfm+i - 
(m -b l)fm = 0 for m > -1, which (using = ifm+i) yields 

3/m + ivlL + + l)/m = 0. (9) 


Comparing with Eq. © allows us to relate the decay exponent a of the velocity amplitude 
to the order m of the Moore-Saffman function. 


a = 


m + 1 
3 


( 10 ) 


Any localized wave motion can be represented as a sum of Moore-Saffman functions of 
different orders m, each leading to a wave component characterized by a distinct decay 
exponent m- This decay exponent agrees with the derivation of Peat^^ for to = 1, and 
with the case m = 0 discussed by Rieutord et aim in the problem of detached layers from 
critical latitudes in a rotating spherical shell. Equation (HUD is also consistent with the 
derivation of Voisinii for internal waves in a stratified fluid (the derivation given in this 
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reference is for the conical wavepacket emitted by a point source, yielding a modified decay 
exponent Gaxi = a + 2/3). 

For a wave beam of order m, we can write explicitly the velocity component along the 
wave beam Ux = 5R({7) as 

= \Uo\ i-j {cm{v)cos{at + a) + Sm,{r])sm{at + a)), (11) 

and the vorticity component ujy = 5(bF) as 

= ^- i-j {-Sm+i{v)cos{at + a) + Cm+i{v)sin{at + a)), (12) 

where t/o = |b^o|e~*“- The argument a accounts for a possible phase shift, through added 
mass effects, between the wave beam oscillation and the source oscillation. These quantities 
CH) and ca are of interest for the experimental measurements based on two-component 
particle imaging velocimetry in the (2f, Z) plane described in Sec. III. 


C. Conservation laws 

We demonstrate now that the order m of the Moore-Saffman function describing a wave 
beam in the far held coincides with the multipole order n of the source from which it is 
emitted. We dehne in the following a source of order n such that the moments of order s of 
the rate of expansion £i{X, Z,t) = V • u are zero for s < n and hnite for s = n. We note hrst 
that the moment of order s of the Moore-Saffman function fm {v) ( 0 ) satishes the property 

(s < m), 

(s = m), (13) 

(s > m), 

so that only the m-th moment of the velocity prohle of order m m is Hnite and non-zero. 
What is needed in addition is to Hnd a conserved quantity involving that moment, and to 
identify m to the multipole order of the source from which the wave beam is emitted. 

Consider Hrst a line monopole (n = 0) along the T-axis, releasing Huid at the Howrate 
q{t) per unit length. The corresponding rate of expansion ^{X^Z,t) has q(t) as its zeroth 
moment, 

q{t) = J £L{X,Z,t)dXdZ, (14) 

and is of the form 

£t{X,Z,t)=qit)S{X)S{Z), (15) 

with 6 the Dirac delta function. For a domain D of boundary C in the (X, Z^)-plane, we 
have, by the divergence theorem. 


V"fmiv)dri = < 


0 

i®7rs! 

oo 


u- ndl 


fi dXdZ, 


ID 


(16) 


where n is the outward normal and dl a positively oriented contour element. We specialize 
to the Hrst quadrant and consider the domain represented in Fig. [21 its boundary starts from 
the origin along a segment Cx of the X-axis, continues with a segment Cz perpendicular to 
the wave beam at a large distance x, and goes back to the origin along a segment Cz of the 
Z-axis. The contributions of Cx and Cz to the contour integral vanish, since the velocity 
is negligible outside the beam. The surface integral is one fourth of the integral over the 
whole plane, owing to the parity of the delta function. We eventually obtain 



Uxdz = \q{t), 


(17) 
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a conservation equation of the type used by Moore and Saffmair^ and Rieutord et al.^ 
involving the zeroth moment of the longitudinal velocity. 

Consider next a line dipole (n = 1), defined such that the flowrate (zeroth moment of the 
rate of expansion /i) is zero, but the first moment 

p{t) = j xn{X,Z,t)dXdZ, (18) 

is finite: it is related to the momentum per unit length imparted to the fluid. The rate of 
expansion /i therefore writes 


fi{X,Z,t) = -p{t)-X{S{X)S{Z)), (19) 

as discussed for example by Piercei^ We write, with Xi an arbitrary coordinate in the plane 
{X, Z) and Ui the associated velocity component. 


V • {xiu) = a;i(V • u) + m. 


( 20 ) 


Integration over an arbitrary domain D yields 


XiU ■ ndl 


XifjL dXdZ - 


: dXdZ, 


ID 


ID 


which becomes, after application to the domain D of Fig. [21 



XiUx dz = \pi + 


Ui dXdZ. 


0 Jo 


( 21 ) 


( 22 ) 


We finally choose Xi = z, the cross-beam coordinate, such that Ui = Uz is negligible every¬ 
where in the quadrant. The last integral vanishes in Eq. (EH) and we obtain 


/ OO 

ZUx dz = jPz{t), 

-OO 


(23) 


a conservation equation involving the first moment of the longitudinal velocity. This equation 
is equivalent to those used by Thomas and Stevenson^S and PeatfiS involving the zeroth 
moments of the pressure and stream function, respectively. We conclude that the velocity 
prohle in a wave beam emitted by a line dipole (n = 1) has a vanishing moment of order 
s = o (HZl and a finite moment of order s = 1 (1231) . 

The previous argument can be generalized to a line source of arbitrary multipole order: the 
wave beam emitted from a source of order n is such that the s-th moment of the longitudinal 
velocity is zero for s < n and finite for s = n. This finite moment writes (see Appendix El) 

/ CO 

z'^Ua:dz = \qz...z{t), (24) 


with qz-.-z the n-th moment of the rate of expansion pL{X, Z,t) along the z axis. Switching 
to complex notation, = 5i(C/oe~*'^*) and qz---z = ^iQz - ze~'‘'^*), and using Eq. (jH]) yield 

J ^"frn[p)dq=\Qz...z. (25) 

Using the property m of fmiv) yields immediately m = n and 


fj — ( Qz-z 

^ nl ’ 


(26) 


conhrming that the order m of the Moore-Saffman wave beam is equal to the multipole 
order n of the source from which it is emitted. 

We can conclude that a beam emitted by a monopole source is essentially an oscillating jet 
of non-zero instantaneous flowrate. On the other hand, a wave beam emitted by a multipolar 
source of order m > 0 contains a set of oscillating shear layers with zero instantaneous 
flowrate, the number of layers in the beam increasing as for large m. The stronger 
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FIG. 3. (Color online) Schematic cross-section of the pulsating (a) and oscillating (b) sources (the 
two extreme states of the oscillation are shown), (c) Vorticity component ojy showing the wave 
beam emitted by the oscillating source, measured by the two cameras with different resolutions 
(fl = 0.68 rad s“^, propagation angle 6 = cos“^((t/2S1) ~ 30°). Different color scales are used for 
the two fields for better visibility. 


shear stress induced by the larger number of layers naturally results in a steeper decay of 
the wave amplitude. 

For a source of arbitrary shape, characterized by an arbitrary multipole expansion, the 
viscous decay of the wave beams in the far-held is dominated by the smallest decay exponent 
a, i.e. by the term of lowest order in the expansion. In practice we can focus on the hrst 
two orders: monopole (m = 0) for any source of Hnite howrate, for which Ux ~ and 

ujy ~ and dipole {m = 1) for a source of zero howrate, for which Ux ~ and 

UJy ^ x~^. 


III. EXPERIMENTAL SETUP AND DATA ANALYSIS 

We have set up an experiment to characterize the inhuence of the multipole order of the 
source on the structure and decay of the inertial wave beam. Measurements are performed 
in a tank of horizontal size Lx x Ly = 150 x 80 cm^, hlled with 50 cm of water, and 
mounted on a 2 m diameter platform rotating around the vertical axis Z. Two wavemakers 
are considered, referred to as pulsating source and oscillating source [see Fig. [31[a,b)]. These 
wavemakers aim to produce effective monopole and dipole sources (properties summarized 
in Tab. [J): 

• The pulsating source consists in a water-hlled horizontal rubber tube, 60 cm long, 
whose volume varies as sin((Tt). The instantaneous radius varies approximately as 
R{t) ~ i?o -b Asin((Tt), with mean radius Rq = 8.9 mm and amplitude A = 0.7 mm 
(the oscillation is harmonic to within A/Rq ~ 8%). This source imposes an oscillating 
how rate per unit length q{t) ~ 27rR qAct cos{at). 

• The oscillating source consists in a horizontal cylinder, 60 cm long, i?o = 3.0 mm in 
radius, whose vertical position Z{t) = Zq + ^sin((Tt) oscillates at frequency a and 
amplitude A = 3.2 mm (Fig. EKb)). This cylinder is a source of zero net howrate 
which can be modeled, at large distances, as an oscillating dipole characterized by 
a dipole moment p(t) = 5i(Pe“*°'‘), with P = [1 -|- C{a)\-KRQAaez and C{a) an 
added mass coefficient. Without background rotation, this coefficient is a real constant 
(7 = 1. With rotation, the coefficient becomes complex and frequency-dependent owing 
to wave generation. The experimental measurement and theoretical determination 
of added mass coefficients have been considered by Ermanyuk and Gavrilov^ and 
Ermanyuk)^ among others, for internal waves. 

The wavemaker (either pulsating or oscillating) is immersed horizontally 10 cm below the 
surface. It is located along the Y axis, at a distance AX = 30 cm from the side wall of the 



















tank. The wavemaker frequency is kept constant, a = 1.18 rad s“^, so that the wavemaker 
velocity Aa is constant. The rotation rate of the platform Q, is varied in the range 0.68 to 
1.68 rad s“^ (6.5 to 16 rpm), resulting in a beam angle 0 = cos“^((7/2n) varying in the 
range 30 — 70°. 

The Reynolds number of the flow in the vicinity of the wavemaker, defined as Re = 
2 RqA<jIv, is Re = 15 for the pulsating source and Re = 23 for the oscillating source. 
Despite the relatively large amplitude ratio (up to A/Rq = 1.06 for the oscillating cylinder), 
these moderate Reynolds numbers indicate that nonlinearities (saturation and generation 
of higher harmonics) can be neglected; see in particular the discussion by Voisin et al^ for 
internal waves, pointing the importance of {A/RQ)Re for saturation. The viscous scale £ ([3]) 
varies in the range 0.6-1.2 mm. For the pulsating cylinder, the ratio Ro/£ ^ 7 —lA indicates 
that the far-field properties of the wave beams are expected at a significant distance xjl 
(the radius of the oscillating cylinder can be made arbitrarily small, but the radius of the 
pulsating cylinder is limited by the design of the rubber tube). 

The two components of the velocity fields (ux,uz) are measured in the vertical plane 
Y = Ly 12 normal to the source axis using a particle image velocimetry (PIV) system 
mounted in the rotating frame. Among the four wave beams emitted by the sources, we 
focus on the one propagating over the longest distance, in the bottom-left direction [see 
Fig-EKc)]- Images of particles are acquired with two 2360 x 1776 pixels cameras operating 
simultaneously with different fields of view. Each PIV acquisition consists in 3 000 image 
pairs (one image per camera) recorded at 3 Hz, which represents 16 fields per source period. 
Cross-correlation between successive images produces velocity fields sampled on a grid of 
295 X 222 vectors, with a spatial resolution of 0.58 mm for the closer view and 2.04 mm for 
the larger view. The combination of the PIV data from the two cameras allows us to resolve 
accurately the spatial scales of the wave field for distances from the source x between 10 mm 
and 1 m. 

Two post-processing steps are applied to the PIV fields. First, the velocity time series are 
phase-averaged at the forcing frequency a in order to filter out contributions from unwanted 
residual flows. These residual flows originate from thermal convection effects (velocities of 
the order of 1 mm s“^), at very small frequencies, and motions due to the precession of the 
rotating platform induced by the Earth rotation^ (~ 0.5 mm s“^), at frequency D. Second, 
we apply a spatial Fourier filter to remove flow structures associated to wave vectors k such 
that kx kz > 0. This procedure is useful to remove secondary wave beams reflecting on the 
tank walls, which intersect the primary beam characterized by kx kz < 0 and induce spatial 
oscillations of the wave envelope [see Fig. He)]. 

Finally, we remap the velocity fields in the tilted frame {x, z) of the beam, with x the 
distance from the source, making the angle 9 = cos“^(cr/2D) to the horizontal, and compute 
the out-of-plane vorticity component ujy. A standard second-order finite difference scheme is 
used to compute uiy, which is comfortably resolved by our twin PIV measurements ensuring 
at least 40 grid points per wavelength at all distance x from the source. The vorticity 
and velocity envelopes of the wave field are finally computed as ujq{x,z) = (2a;^)^/^ and 
uo{x,z) = (2u^)^/^, with 0 a temporal average. 


IV. RESULTS 

We first compare the spatial decay of the vorticity envelope of the wave beam for the 
pulsating and oscillating sources. Vorticity is used here instead of velocity to compare 
against the theory m because it is less sensitive to residual large scale flows. The centerline 
vorticity ujo{x, z = 0) is plotted as a function of the distance x from the source in Fig. Ha') foi' 


Source 

Ro 

[mm] 

A 

[mm] 

Aa 

[mm s“'^] 

Re 

Vorticity dec 
Theory 

ay exponent 
Experiment 

Pulsating (m = 0) 

8.9 

0.7 

0.82 

15 

-2/3 

-0.64 ±0.09 

Oscillating (m = 1) 

3.0 

3.2 

3.77 

23 

-1 

-0.99 ±0.05 


TABLE I. Properties of the pulsating and oscillating sources: mean radius J?o, oscillation amplitude 
A, velocity amplitude Aa, Reynolds number Re = ^R^Aajv (with v — 10~® re? s“^ the kinematic 
viscosity of water). The vorticity decay exponents correspond to the average (± standard deviation) 
over the four experiments at different cr/2D [see Fig. Hb)]- 
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FIG. 4. (Color online) (a) Vorticity envelope ojo{x,z = 0) along the direction of the wave beam 
for the pulsating (thick line) and oscillating (thin line) sources (fl = 1.05 rad s“^, 6 ~ 56°, £ = 
0.76 mm). The lines show power laws of exponents —2/3 and —1 expected theoretically for a 
monopole and a dipole source respectively, (b) Vorticity decay exponent as a function of the 
propagation angle 9 = cos~^ ((T/2fl): (o), pulsating source; (□), oscillating source. 



n 

FIG. 5. (Color online) Normalized transverse vorticity profiles ujy{x,r],t)/Loo{x,0) as a function of 
the reduced transverse coordinate r] = compared to the theoretical vorticity profiles 

Cm+i(7?)/cm+i(0). O, pulsatiug source (m = 0); □, oscillating source (m = 1). The profiles are 
measured at a distance xj£ = 330 from the source, in the region where both beams show a power 
law decay, and are shown for the time at which the vorticity is maximum at the center of the beam. 


n = 1.05 rad iO ~ 56°). The distance is normalized by the viscous length £, and vorticity 
Wo is normalized by the source velocity Aa and i. For both sources a power law decay of the 
vorticity is observed far from the source, but with a different exponent: wq ~ ^.-i.ooio.os 
for the oscillating source, starting close to the source [xjl > 10) and extending over nearly 
two decades; wq ~ pulsating source, starting much further from the 

source {x/£ > 200) and hence visible over a limited range of x. These exponents are in 
good agreement with the predictions for a monopole and a dipole source, x~'^l^ and x~^ 
respectively [see Eq. (inD]. Exponents measured at other values of cr/2fl = cos0, reported 
in Fig.|4Kb), are consistent with these numbers and do not show any trend with 0. 

Another distinctive property of a wave beam emitted by a monopole or a dipole distur¬ 
bance is the transverse vorticity profile, given by the function fi(r/) or / 2 (? 7 ) respectively. 
Figure |S] shows the normalized vorticity profile ujy(x, z, t)/uJo(x, z = 0) for both sources as a 
function of the reduced transverse coordinate r/ = at the time at which ujy is 

maximum at the centerline. The normalized profiles, shown here for a distance xjl = 330, 
are independent of x provided that x is large enough, i.e. in the region showing a well-defined 
power-law decay {xj I > 200 for the pulsating source and xj I > 10 for the oscillating source). 
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FIG. 6. (Color online) Contour plot of the vorticity envelope ujq from the close-view camera for 
the pulsating source at 9 = 56°, showing the transition from a bimodal to a unimodal wave beam. 
The white circle shows the source and the white line the far-held wave beam axis. The star marker 
indicates the location Xcji = 110 at which t<Jo(a;) is maximum in Fig. |4](a). 


We compare these profiles against the Moore-Saffman functions taken at the same phase, 
ujy{x,r],t)/ujQ{x,r] = 0) = Cm-i-i(? 7 )/cm+i( 0 ), for m = 0 and m = 1. The agreement is excel¬ 
lent for both sources, for \ri\ up to 5, clearly confirming that the order of the Moore-Saffman 
function that best describes the wave envelope is governed by the multipole order of the 
source. At larger distance from the beam centerline {\r]\ > 5), the discrepancy between the 
theoretical and experimental profiles probably originates from residual fluid motions associ¬ 
ated to reflected wave beams that cannot be eliminated by the Fourier filtering procedure. 

The flowrate across the wave beam provides another confirmation of the match between 
the wave beam order and the source multipole order. We compute the flowrate q{t) per 
unit length by integration over rj of the instantaneous velocity Ux{x,r]), using a truncation 
at \r]\ = 6 to reduce disturbance from fluid motions out of the primary beam of interest. 
The theoretical flowrate amplitude for each wave beam is qth = 77 = 0) 

for the monopole source and zero for the dipole source. We find a normalized flowrate 
q/qth — 1-0 ± 0.05 for the pulsating cylinder, and 0.07 ± 0.05 for the oscillating cylinder 
(the non-zero value in the latter case originates from the truncation of the integral). This 
normalized flowrate does not depend significantly on the distance from the source in the 
far-held for both disturbances (to within ±0.05). 

We Hnally turn to the description of the distance x beyond which the scaling law holds for 
the far-held decay of the wave amplitude. Figure SKa) shows that the vorticity amplitude 
decreases at all x for the oscillating source, with a well-dehned power law beyond x(£ 10 , 

while it is non-monotonic for the monopole source with a maximum at xd£ ~ 110 ± 10 . 
This non-monotonic prohle originates from the transition from a bimodal beam close to the 
source to a unimodal beam far from the source, as illustrated by the close-up view in Fig. [5] 
Because of the large extent of the pulsating cylinder compared to the viscous length, the 
wave held close to the source corresponds to two separate beams propagating in the same 
direction 0 = cos“^(( 7 / 2 fl), with a shift of z = ±i?o = ±8.9 mm from the axis of the far-held 
beam. Close to the source, the transverse prohle is therefore bimodal, with a local minimum 
at z = 0 , while it becomes unimodal for large distances with a maximum at z = 0 , with a 
transition occurring around Xc- 

This bimodal-to-unimodal transition, a classical feature of waves emitted from a source 
of large extent, has been mainly described for internal waves excited by oscillating 
disturbancesFigure [IJa,b) provides a simple sketch of this transition for a pul¬ 
sating and an oscillating cylinder of large extent. The oscillating boundary layer over the 
cylinder detaches at the four critical lines, where the local slope equals the wave beam angle 
0, forming eight wave beams of non-zero instantaneous howrate. The far-held beam in a given 
direction results from the merging of two parent beams propagating along the same direction. 
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that are in phase for the pulsating cylinder and out-of-phase for the oscillating cylinder. The 
combination of two in-phase parent beams of order m separated by a normalized transverse 
distance rjo is given by fm{v + ? 7 o/ 2 ) + fniiv — ^o/2), which far from the source {rjo <C 1) 
simply gives a beam of order m. On the other hand, the combination of two out-of-phase 
parent beams of order m is given by /miv + ? 7 o/ 2 ) - fm{r] - rjo/2) ~ mf'miv) = *%/m-i-i(b). 
yielding a beam of order m -I- 1: this is consistent with the observation of a far-field wave 
beam of zero flowrate (m = 1) emitted from the two nonzero flowrate sources {m = 0) at 
the critical lines for the oscillating cylinder. 

Determining the merging distance Xc requires the full resolution of the flow close to the 
wavemaker. This computation is given in Hurley and Keady^ and Voisin^ for oscillat¬ 
ing cylinders and spheres in a stratified fluid. Qualitatively, we can estimate Xc from the 
spreading law 6{x) ^ of each parent beam, yielding a merging distance Xc such that 

(^(xc) — 2i?o given by x^ji ~ ■ The radius ratio between the pulsating and oscillating 

cylinders {Rq = 8.9 and 3.0 mm, respectively) indicates that xd^ is expected much larger 
for the pulsating cylinder. For the oscillating cylinder the predicted merging distance is of 
order of the cylinder radius, so that the bimodal wave beam cannot be observed, which is 
consistent with the monotonic decay of ojo(x). 


V. CONCLUSION 

In this paper we analyzed the decay of a two-dimensional inertial wave beam and showed 
that it is set by the multipole order n of the source from which it is emitted. Experimental 
measurements are reported for sources of the first two orders: monopole (n = 0) and dipole 
(n = 1). We find that the structure of the wave beam is well represented by a Moore- 
Saffman^i (or Thomas- Stevensort^) function of order equal to the multipole order of the 
source. The wave envelope decays as a power law of the distance from the source, with 
an exponent governed by the order of the source for the velocity and 

for the vorticity). These properties, demonstrated here for inertial waves in rotating fluids, 
should also hold for internal waves in stratified fluids. 

The steeper decay of the wave amplitude as the multipole order is increased indicates that 
the far-held structure of a wave beam is dominated by the lowest order of the source, i.e. by 
its Hrst nonzero moment (flowrate, momentum or higher order moment). In most practical 
situations, the nature of the source can be discriminated by its instantaneous howrate, either 
nonzero for a monopole source or zero for a dipole or higher order source. This means that 
a wave beam emitted from any source of nonzero instantaneous flowrate must be dominated 
in the far held by its monopole component. Such monopole sources are relevant to most 
natural hows (e.g. conical wave beams in the huid core of planets, internal tide generation 
over ocean topography): the detachment of the oscillating boundary layers at the critical 
lines produces oscillating jets in the bulk, corresponding to a weak spatial decay. On the 
other hand, oscillating disturbances immersed in the huid (a typical conhguration of most 
laboratory experiments in rotating or stratihed huids) produce in the far held beams of zero 
instantaneous howrate, resulting in a stronger spatial decay. 
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Appendix A: Conservation laws for arbitrary multipole order 

Equations dm and (1251) relate the zeroth and hrst moments of the longitudinal velocity 
to the howrate and momentum of monopole (n = 0) and dipole (n = I) sources. In 
this appendix we generalize these relations to sources of arbitrary multipole order n. We 
introduce arbitrary orthogonal coordinates {xi,X 2 ) in the plane {X,Z), with (ui,U 2 ) the 
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associated velocity components. The rate of expansion fi{X, Z, t) of a source of multipole 
order n has zero moments of order 0 to n — 1 and finite n-th moments given by 

lii-iAi) = j Xi^ ■ ■ ■ Xi^fi{X,Z,t)dXdZ. (Al) 

These n-th moments are tensors of rank n composed of n -|- 1 independent scalars. The rate 
of expansion is of the form 


I'—l'l" 2 2 „„ 

(A2) 

*1 = 1 *n=l 

Using the identity 

V • (xji • • -Xi^u) • u) -I- Ui^Xi^ • • • H - Vxi^--- Xi^_^Ui^, (A3) 

and integrating over the domain D of Fig. [2] yield 



■Xi^Uxdz = 


n oo 

{ui^Xi^ - \-Xi^ •••Xi^_-,Ui^)dXdZ. 


(A4) 


We choose {xi^X 2 ) = (x^z) and zi = • • • = = 2, so that Xi^ = ..• = = z. The 

corresponding transverse velocity Ui^ = • • • = = Uz is negligible, so the last integral 

vanishes in Eq. (|A4p . We finally obtain Eq. (P4|l . which generalizes Eqs. (fT7|l and (P5|) to 

arbitrary multipole order n. 
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